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ABSTRACT 

We present Keck LRIS spectroscopy of the black hole-hosting globular cluster RZ2109 in the Virgo ellip- 
tical galaxy NGC 4472. We find that this object has extraordinarily broad [O III]A5007 and [O III]A4959 
emission lines, with velocity widths of approximately 2,000kms~^. This result has significant implica- 
tions for the nature of this accreting black-hole system and the mass of the globular cluster black hole. 
We show that the broad [O III]A5007 emission must arise from material driven at high velocity from the 
black hole system. This is because the volume available near the black hole is too small by many orders 
of magnitude to have enough [OIII] emitting atoms to account for the observed L([0 III]A5007) at high 
velocities, even if this volume is filled with Oxygen at the critical density for [O III]A5007. The Balmer 
emission is also weak, indicating the observed [OIII] is not due to shocks. We therefore conclude that 
the [O III]AA4959,5007 is produced by photoionization of material driven across the cluster. The only 
known way to drive significant material at high velocity is for a system accreting mass near or above its 
Eddington limit, which indicates a stellar mass black hole. Since it is dynamically implausible to form 
an accreting stellar mass black hole system in a globular cluster with an intermediate mass black hole 
(IMBH), it appears this massive globular cluster does not have an IMBH. We discuss further tests of 
this conclusion, and its implications for the Mbh — -^stcUar and Mbh ~ cr relations. 
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1. INTRODUCTION 

Maccarone et al. (2007) presented the first unambigu- 
ous evidence for a black hole in a globular cluster. This 
was based on an XMM observation of strong variability in 
a high luminosity X-ray source (~ 4 x 10'^^ ergs/s) located 
in the spectroscopically confirmed globular cluster RZ2109 
in the Virgo elliptical galaxy NGC 4472. The high X-ray 
luminosity, about an order of magnitude greater than the 
Eddington luminosity of a neutron star, requires either a 
black hole or multiple neutron stars in the old stellar pop- 
ulation of the globular cluster. The strong variability rules 
out multiple neutron stars as the source of the bright X-ray 
emission, and thus indicates the globular cluster RZ2019 
hosts a black hole X-ray source. 

Subsequently, Zepf et al. (2007, hereafter Z2007) dis- 
covered broad [O III]A5007 emission from this black hole- 
hosting globular cluster. Specifically, Z2007 analyzed ex- 
isting VLT spectra of this object covering a wavelength 
range of 5000 — 5800 A. They identified a broad emission 
line in these data at a wavelength in very close agreement 
(within 10 kms"^) with that expected for [O III]A5007 at 
the radial velocity of the globular cluster previously deter- 
mined from its stellar absorption lines. Z2007 also found 



that the [O III]A5007 emission line was much broader than 
any possible velocity dispersion for the globular cluster, 
with an estimated width of several hundred kms~^, and 
the possibility of much broader wings in the data. They 
noted such a large velocity width clearly points to the black 
hole as the source of the power driving the nebular emis- 
sion. 

As described in Z2007, understanding how the broad 
[O III]A5007 emission is produced provides a way to con- 
strain the mass of the black hole in this globular cluster. 
One possibility is that the black hole is a stellar mass sys- 
tem, accreting near or above its Eddington limit. Such 
systems can drive strong, high velocity winds (Begelman, 
King, & Pringle 2006 and references therein). In this 
case, broad [OIII] could be produced either by shocks as 
the wind passes through the globular cluster's interstel- 
lar medium, or by photoionization of this wind material 
across the cluster. Alternatively, it could be an interme- 
diate mass black hole, in which case it would be far below 
its Eddington luminosity given the observed X-ray lumi- 
nosity. In this case, no strong winds would be driven, and 
any broad [OIII] emission would have to be produced by 
photoionizing of material in very close proximity to the 
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black hole. 

The X-ray data alone do not clearly distinguish between 
these possibilities. There are also a wide range of theoret- 
ical predictions for the fate and properties of black holes 
in globular clusters. For example, extrapolations of the 
M^BH — o and the Mbh — Mgteiiar relations tend to predict 
IMBHs in globular clusters (see Section 4), and some mod- 
els find sufficiently dense young globular clusters cores may 
experience runaway stellar mergers that may form IMBHs 
(e.g. Portegies Zwart & McMillan 2002, Frcitag ct al. 2006, 
but see also Yungelson et al. 2008). On the other hand, 
as stellar mass black holes sink to the centers of globular 
clusters through mass segregation, dynamical interactions 
between them may eject most or all of the black holes from 
the core, cither into the outskirts of the globular cluster 
(Mackey et al. 2007), or out of the cluster completely (e.g. 
Sigurddson & Hernquist 1993, Kulkarni et al. 1993). 

In this Letter, we present new optical spectroscopy 
aimed at understanding the physical nature of the accret- 
ing black hole source in the globular cluster RZ2109. The 
optical spectra used in Z2007 were available because of an 
ongoing kinematic study of the NGC 4472 globular clus- 
ter system using FLAMES, a multi- fiber spectrograph on 
the VLT. The FLAMES spectra have very limited wave- 
length coverage and a modest signal-to- noise ratio. The 
work presented here both significantly increases the wave- 
length coverage to investigate diagnostic line ratios, and 
improves the signal to noise ratio to better constrain the 
width of the broad [O III]A5007 emission line. 

2. OBSERVATIONS 

We obtained optical spectra of the NGC 4472 globular 
cluster RZ2019, which contains the X-ray source XMMU 
1229397+07533, on UT 2007 December 17-18 with the 
Low Resolution Imaging Spectrometer (LRIS; Oke et al. 
1995) on the Keck I telescope. The observations, taken at 
a position angle of 74.4 deg, through a 1'.'5 slit, used the 
D560 dichroic, the 400 lines mm~^ grism (Abiazc = 3400A) 
on the blue side, and the 400 lines mm^^ grating (Abiazc = 
8500A) on the red side. The spectral resolution for objects 
filling the slit was R ~ 400 and 700, on the blue and red 
sides respectively. The nights were not photometric, and 
the seeing ranged from moderate (1'.'2; first night) to poor 
(1'.'6; second night). The exposure times were 1500 s each 
night. 

Data processing followed standard slit spectroscopy pro- 
cedures. In particular, we extracted the spectra using a 
1'.'5 wide aperture and calculated the spectral dispersions 
using lamp spectra obtained earlier in the nights, adjust- 
ing the wavelength zero point based on telluric emission 
lines. The spectra were flux calibrated using observations 
of standard stars from Massey & Gronwall (1990) obtained 
during the same nights. Since the nights were not photo- 
metric, the spectra were scaled to match the observed op- 
tical brightness of RZ2109 (F = 21.0). We also corrected 
the atmospheric A and B band bands based on the stan- 
dard star spectra. The final, calibrated spectrum is shown 
in Figure 1. 

3. RESULTS 

In this letter, we focus on two key aspects of the re- 
markable spectrum of RZ2109 shown in Figure 1. These 



are - 1) the very broad velocity width of the [O III]AA4959, 
5007, and 2) the [O III]A5007 line luminosity and the large 
ratio of the [O III]A5007 to H/3. We also note that this 
spectrum clearly shows [O III]A5007, [O III]A4959, and 
likely [O III]A4363 in emission at the radial velocity of the 
globular cluster previously determined from its stellar ab- 
sorption lines. This confirms beyond any doubt that these 
emission lines are from the globular cluster RZ2109, as 
found by Z2007 from the [O III]A5007 line alone. 

3.1. The Extraordinary [O III]XX4959, 5007 Velocity 

Width 

The spectrum in Figure 1 shows that the [O III]A4959 
and A5007 emission lines are extremely broad, so broad 
that the lines are blended. To determine the velocity width 
of these lines we first use standard deblending procedures. 
Specifically, we fit the two lines with two Gaussians, keep- 
ing the width of the lines the same because they come 
from the same initial atomic state. This gives a FWHM 
of the lines of 33 A, corresponding to a velocity FWHM of 
approximately 2000 kms^^. 

Figure 1 shows the lines have stronger central peaks 
and broader wings than the best fitting Gaussian func- 
tions. We therefore also calculate the velocity width of the 
lines by determining the longest wavelength at which the 
[O III]A5007 line clearly has significant flux. We determine 
the corresponding shortest wavelength for the [O III]A4959 
line to estimate the lowest velocity component. Using this 
technique, we find a velocity extent of the emission of 
at least +1600 kms"^ and -1400 kms~^ Whether the 
[O III]AA4959, 5007 lines are described as having FWHM 
of 33 A, or having a velocity extent of at least ±1500 
kms~^, these lines are extraordinarily broad, as broad as 
almost any [O III]A5007 lines observed, while being located 
in a globular cluster with an estimated mass of a few mil- 
lion solar masses and an escape velocity several orders of 
magnitude less than the observed velocity width. 

3.2. Line Luminosities and Ratios 

The [O III]A5007 fine luminosity and its comparison to 
other diagnostic lines such as H/3 place important con- 
straints on the nature of the accreting black hole system. 
One way to calculate L([0 III]A5007) is to use the Gaussian 
fits to the [OIII] lines described above. These give fluxes of 
1.8 X 10^^^ ergss~-^cm~^ and 4.4 x 10"-^^ ergss~-^cm^^ for 
[OIII]4949 and [O III]A5007 respectively. Adopting a 16 
Mpc distance for NGC 4472, then gives a line luminosity 
of L([0 III]A5007)= 1.4 x lO^^ ergss"!. 

An alternative approach to determining the [O III] A5007 
luminosity that avoids Gaussian fitting is to determine 
the total [OIII]A4959 + A5007 luminosity, and adopt the 3 
to 1 ratio of L([0 III]A5007) to L([0 III]A4959) given by 
atomic physics (e.g. Dimitrijevic et al. 2006, Storey & Zeip- 
pen 2000). We do this by measuring the straightforward 
equivalent width and total [0III]A4959 + A5007 flux, and 
then using the above ratio between the lines to determine 
L([0 III]A5007). This calculation gives an [O III]A5007 
equivalent width of 31 A, and L([0 III]A5007) = 1.4 x lO^'^ 
ergss^^. While the close agreement of this determination 
with that above is undoubtedly partially fortuitous, it is 
also reassuring that different techniques for determining 
the line luminosities give similar answers. We also note 
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that the L([0 III]A5007) found here is significantly higher 
than the best estimate in Z2007. Z2007 noted that the Une 
might have broad wings, and that if the broad wings were 
real, their velocity width and L([0 III]A5007) estimates 
would be lower limits. 

The [O III]A5007 to H/3 ratio is a valuable diagnostic 
emission line mechanism. It is immediately apparent from 
Figure 1 that the Balmer line emission is not strong; the 
Balmer lines are all seen in absorption. Balmer absorp- 
tion lines are expected from the stellar population of the 
globular cluster. The H/3 emission can be estimated by 
comparing the amount of H/3 absorption expected from 
stellar population models to that observed, attributing 
any difference to filling in of the stellar absorption line 
by H/3 emission from the black hole system. To carry out 
this measurement, we use the stellar populations models 
of Bruzual & Chariot (2003), and those of Vazdekis et al. 
(2008) based on the spectral library of Sanchez-Blazquez 
et al. (2006), since these provide spectra at fairly high 
resolution. For each set of models, we adopt an age of 12 
Gyr and set the model metallicity to match the observed 
optical colors of RZ2019. We then analyze the spectral 
output of the models in the same way as the data. Both 
models predict an equivalent width of 2.6 A for the H/3 ab- 
sorption line, while the same measurement of the observed 
spectrum gives 1.5 A. 

Our analysis therefore suggests an H/3 equivalent width 
of about 1.1 A giving a very rough L[0 IIl]A5007/L[H,g] 
ratio of 30. We caution that this is somewhat dependent on 
the stellar populations modeling of the stellar Balmer ab- 
sorption. There are also hints in the current data that the 
Balmer emission may not be as broad as the [O III]AA4959, 
5007 lines, and we plan to address this and other questions 
about the detailed shape of the [O III]AA4959, 5007 line 
complex in future work with higher spectral resolution, we 
will address this further in a future paper. Here wc simply 
note that the [O III]A5007/H/3 ratio appears to be large. 

4. DISCUSSION 

Zepf et al. (2007) noted three possibilities for produc- 
ing the broad [O III]A5007 emission. One is that it could 
be produced by shocks from a strong, high velocity wind 
driven from a black hole. Driving such winds requires 
a system undergoing mass transfer at or above its Ed- 
dington limit, which given the observed Lx ~ 4 x 10^^ 
ergss"^, would indicate a stellar mass black hole system 
of Mbh — 10 M0. Alternatively, the velocity width could 
be generated by a strong wind driven across the globular 
cluster, but the ionization source could be photoionization 
from the powerful central system. Because this case in- 
volves strong winds, it too implicates a stellar mass black 
hole system. The third possibility is that the high ve- 
locities are due to gravitational motions very close to the 
central black hole, and the material is photoionized. This 
scenario likely implicates an IMBH to more readily give 
the high velocities at sufficient distances to allow [OIII] 
emission. 

The new resiilts presented here - the very broad 
[O III]A4959 and A5007 emission lines extending over more 
than one thousand kms~^, and the large [O III]A5007/H/3 
ratio - strongly constrain these possibilities. First, the 
large [OIII] /H/3 ratio indicates the [O III]A5007 is not pro- 



duced directly by shocks. Secondly, the large line widths 
extending to high velocities rule out the idea the velocity 
widths seen in the [O III]AA4959, 5007 lines can be due to 
gravitational motions very near a black hole. This is be- 
caiise the vohime available near the black hole where the 
line widths would be as large as observed is too small by 
orders of magnitude to have enough [OIII] emitting atoms. 

To calculate the maximum possible [O III]A5007 flux for 
material orbiting around a black hole, we first find the 
volume available in a shell corresponding to a given cir- 
cular velocity around the black hole. We then determine 
the maximum possible [O III]A5007 flux at this velocity 
by assuming this volume is filled completely with mate- 
rial at the critical density of [O III]A5007 with a solar 
O/H ratio, and taking the coUisional excitation rate for 
[O III]A5007 (Osterbrock & Ferland 2006). We find that 
the maximum flux for material with a velocity of 1500 
kms-^around a 1000 M© black hole is Fx{[0 III]A5007) 
= 8 X 10~^^ ergss"^ cm~^A~^. This is the maximal flux 
in that it assumes the entire available volume is at the 
critical density of [O III]A5007, that all the O atoms are 
doubly ionized, and that the O/H ratio is solar, while 
the best estimate for the metallicitj^ of RZ2109 is about 
1/50 Zq (Maccarone et al. 2007, Rhode & Zepf 2001). 
In contrast to this maximum flux from material orbiting 
around a black hole, the observed ^a([0 III]A5007) at ve- 
locities around 1500 kms~^ is somewhat larger than 10~^* 
ergss"^ cm~^A~^ (see Figure 1), or more than 10^ times 
greater. Thus, even the maximal emission from material 
orbiting around a black hole fails by many orders of magni- 
tude to match the flux in the observed broad [O III]A5007 
line. One can increase the flux by assuming an even higher 
O abundance, or by increasing the available volume by 
going to higher Mbh or lower velocity. However, the dif- 
ference is so great and the presence of signiflcant flux at 
large velocities clearly observed, that these do not begin 
to reconcile this scenario with the observations. 

Therefore, the only scenario consistent with the data is 
that in which the black hole system drives material across 
the cluster, and that this material is also photoionized by 
the central X-ray source. Such a scenario strongly im- 
plies that the black hole is a stellar mass black hole and 
not an IMBH, based on both theoretical and observational 
considerations. Extensive theoretical work has found that 
strong winds like those found here are driven only by sys- 
tems with mass accretion rates similar to or larger than the 
Eddington limit (e.g. Proga 2007 and references therein). 
Given the observed X-ray luminosity of 4 x 10^^ ergss"^, 
this strongly implies a stellar mass black hole, with a mass 
of Mbh — 10 giving a good fit to the X-ray data in 
this scenario (Maccarone et al. 2007). 

A stellar mass for the accreting black hole system in 
RZ2109 is also supported by the high L/LEdd found for 
the small fraction of cxtragalactic AGN with observed 
broad velocity widths in [O III]A5007. One well-studied 
example is NGC 1068, which has long been known to have 
highly blueshifted [O III]A5007 components with the emis- 
sion dominated by photoionization (Groves et al. 2004, 
Kinkhabwala et al. 2002), and for which the central source 
appears to be above its formal Eddington limit (e.g. Bock 
et al. 2000). A small fraction of narrow- line Seyfert 1 
galaxies are also observed to have either an underlying 
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broad or significantly blueshifted [O III]A5007 component, 
which is found to be related to high i/LEdd ratios (e.g. 
Komossa et al. 2008, Boroson 2005), further suggesting 
RZ2109 is near its Eddington luminosity. 

Thus, theory and observation indicate that the black 
hole in this globular cluster system is very likely a stellar 
mass black hole. As discussed in Z2007, it seems highly 
unlikely that a close, accreting black hole binary system 
could form in a globular cluster with an IMBH at its cen- 
ter. This is because the dynamical processes by which 
close black hole X-ray binaries form in globular clusters 
happen in their cores, but dynamically, a globular cluster 
core containing an IMBH will have typically ejected its 
stellar mass black holes (see Z2007 for full discussion). 

The globular cluster RZ2109 appears to have only stel- 
lar mass black hole(s), with no intermediate mass black 
hole. As such it falls far off the relations between black 
hole mass and the mass and velocity dispersion of the stel- 
lar system in which it is found. To show this explicitly, 
we first calculate the Mbh for RZ2109 implied by the 
Mbh - AfstcUar relation found by Haring & Rix (2004). 
For the stellar mass of RZ2109, we use its absolute lu- 
minosity of My — —10, and a {M/L)v = 2 standard 
for globular clusters to estimate a stellar mass of about 
2 X 10^ Mq. Based on the Haring & Rix (2004) relation, 
a Mbh of 900 Af© would be predicted for RZ2109, with 
an uncertainty of a factor of two. This is far greater than 
the stellar mass of ^ 10 M© indicated by the calculations 
above. Extrapolating the Mbh — f relation to RZ2109 
tends to predict even larger black hole masses than the 
Mbh — AfstoUar relation. A direct determination of the a 
of RZ2109 is not feasible, but we can compare it to the 
a found for Galactic globular clusters of similar absolute 
magnitude. These are the luminous clusters such as 47 
Tuc, Omega Cen, and M54, which have a ranging from 
about 12 — 20 kms^^. Adopting a typical value of 16 



kms~^, the Tremaine et al. (2002) AIbh ~ c relation gives 
a predicted central black hole mass in RZ2109 of 5000 Mq. 

The present analysis would seem to indicate that not 
all old, metal-poor stellar systems form black holes con- 
sistent with Mbh — Afsteiiar and Mbh — c relations found 
for more massive galaxies. It also indicates that stellar 
dynamics did not lead to coalesence into an intermediate 
mass black hole, at least in this one fairly massive glob- 
ular cluster. More speculatively, the results suggest that 
the dark matter halos that presumably surround the early 
galaxies which grow hierarchically to make massive galax- 
ies at the current epoch are important for retaining gas 
which is subsequently accreted to fuel the BH growth and 
set the Mbh — cr relation. To further test these conclu- 
sions, additional work on forming accreting stellar mass 
black hole systems in globular clusters with and without 
an IMBH would be valuable, as would further understand- 
ing of winds driven from systems undergoing mass trans- 
fer near or above their Eddington limit. Data with higher 
spectral resolution would allow a more detailed study of 
the velocity structure in the [O III]AA4959, 5007 emission 
lines which would help constrain the geometry of the emit- 
ting material. We also note that the extended nature of 
the [O III]A5007 emission in RZ2109 can be directly tested 
by spatially resolved spectroscopy with HST. 
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Fig. 1.— A plot of Fx vs. A for RZ2109, the NGC 4472 globular cluster in which the black hole system XMMU 1229397+075333 is found. 
This figure shows the remarkably strong and broad [O IIIJAA4959, 5007 [OIII] emission lines. The upper inset figure shows the [O IIIJAA4959, 
5007 region in more detail. The light dashed lines on the inset figure arc the Gaussian fits to the [O IIIJAA4959, 5007 lines, with a FWHM 
of 2,000 knis^-'. This demonstrates the extraordinary width of these lines. The paper focuses on these lines and the implication of their 
luminosity and velocity width. There also appears to be velocity structure beyond the individual Gaussian fits, but higher spectral resolution 
than the R ~ 400 of these data is needed for a reliable assessment of structure within the very broad lines. The gap around 5600 A is due 
to the dichroic in the LRIS spectrograph. The stellar absorption lines are from the stellar population of the host globular cluster, with some 
likely filling in of the Balmer lines due to weak Balmer emission, as discussed in the text. 



